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Abstract 



We consider the wave equation with a damping term on a partially rectangular planar do- 
On , main, assuming that the damping is concentrated close to the non-rectangular part of the domain. 
■ Polynomial decay estimates for the energy of the solution are established. 
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; 1 Introduction and statement of result 

7-H ■ 

The purpose of this note is to show how the methods of [I], [3] apply to estimate the energy decay 
rates for the damped wave equation on a class of planar domains, including some ergodic billiards. 
In a situation when the geometric control condition of Bardos, Lebeau, and Rauch does not hold, 
we obtain a polynomial decay estimate for the energy of the damped wave, with respect to a stronger 
norm of the initial data. In order to formulate the main result, we shall begin by recalling some 
standard notation and assumptions. 

Let Q, C R n , n > 2, be a bounded connected domain, with dQ, G C°°. When < a G L°°(f2) is a 
non-negative smooth function on 17, we consider the following initial-boundary value problem, 



(-Z) t 2 - A + 2ia(x)D t ) u = 0, (t, x) G R+ x U, 

u = on R+ x dn, (1.1) 
u|t=o = «o£ H^(n), D t u\ t =o = ux G L 2 (n). 

Here Dt = d%ji and we shall assume throughout that the damping coefficient a does not vanish 
identically. 

Associated with the evolution problem l|l.l|) is the solution operator U(t) = e lt-A , t > 0, acting in the 
Hilbert space of the Cauchy data TL = Hq{Q) x L 2 (Q) and mapping (uq,Ui) G TL to (u(t, •), D t u(t, •)). 
Here we equip 7i with the norm 



"o 



Vn His + || ni 1 1^2. 
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The infinitesimal generator A of the semigroup U(t) is the operator 



A=( °. ,, * ) (1.2) 
y —A lia[x) J 

with the domain D {A) = (H 2 n ifg) x //q. It follows that the spectrum of A is discrete, and from jS] 
we recall that if A G C is an eigenvalue of A then < Im A < 2| | a \ . 

For future reference, let us notice that when U G D(A), then Im (AU, U) n > and hence 

II (A-^r 1 \\c(H,H) ^J^xy ImA<0. (1.3) 

In this note we shall be concerned with the energy of the solution u(x, t) of (|1-1|) at time t, 

E(u, t) = \J (l V^| 2 + | A^ 2 ) dx. (1.4) 

It is easily seen that E(u, t) is nonincreasing as t — > oo and from [H] we may also recall that E(u, t) — ► 
as £ — ► oo, for each (uo,tti) G Under the geometric control condition stating that there exists 
To > such any billiard trajectory in £1 of length > To meets the open set {x;a(x) > 0}, the uniform 
exponential decay of the energy has been established by Bardos, Lebeau, and Rauch pQ. In the general 
case, without any assumptions on the underlying dynamics, it has been proved by Lebeau 8 that the 
decay rate of the energy is always logarithmic, provided that the initial data in (jl.l() are measured 
with respect to a stronger norm. In this note, we shall derive a polynomial decay estimate for the 
energy, for a class of planar domains and damping regions, in the case when the geometric control 
condition of P fails to hold. 

We shall consider the class of partially rectangular domains £1 C R 2 . By this we mean that 0, is 
connected, has a boundary that is piecewise C°°, and contains a rectangle R C f2, such that if we 
decompose the boundary of R into pairs of parallel segments, dR = Ti Ul^, then Tj C d£l for at least 
one j, say for j = 1. We shall write f2 = R U W, where W is the non-rectangular part of ft. 

Example. The Bunimovich stadium S, defined as the union of a rectangle R = {(x,y);x G [0,1], y G 
[0,2/3]}, (3 > 0, with two semicircular regions centered at (0, (3) and (1,(3), with radius (3, which lie 
outside of R, is a partially rectangular domain. It is well-known [2] that the geodesic flow in S, obeying 
the law of reflection at the boundary, is ergodic. In what follows, the ergodicity of the underlying 
classical flow in O will not play any special role in our considerations. 

Theorem 1.1 Let Q, C R 2 be a partially rectangular domain, O = R U W, and let < a G C(Q) be 
such that a > in W , the closure of W in £1. Then for each k > there exists a constant > 
such that for each (uq,ui) G D(A k ) we have 



E{u,tf/ 2 < C fc (1 ° gt ^ ■ || (u ,ux) \\ D(Ak) , t>2. (1.5) 

£2 



+1 



Remark. We notice that in the case when the damping coefficient a vanishes outside of a neighborhood 
of W in $7, the geometric control condition of fails, due to the existence of the invariant set for the 
classical flow constituted by the bouncing ball orbits parallel to the pair of segments T2 and staying 
within the rectangle R. Therefore (see pQ), no uniform decay estimate for the energy with respect to 
the energy norm of the initial data is possible. 
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Theorem 1.1 is proved in section 2. We shall finish this section by briefly discussing the issue of 
optimality of the bound (|1.5|) . First, arguing heuristically on the classical level, let us notice that 
every ray starting in the interior of the rectangle and forming an angle of size ^, i > 1, with the 
direction of will reach the non-rectangular part of the domain after a time which is of the order ~ k. 
When discussing the quantum picture, for simplicity, we shall replace the rectangle R by an infinite 
strip of the form R x x [0, ir] y and consider the wave evolution of the functions e^{x, y) = p(x) sin ky, 
p G <S(R), tp £ Cq°(R), II p ||^2 =1. Here (p(£) = f e~ tx ^p{x) dx is the Fourier transform of p. Then 
e/c is an 0(l)-quasimode for the Dirichlet Laplacian Ad in H x x [0, Tv] y , and a simple calculation shows 
that 

^cos t a/— A~d^ efc = (cos tk) + O , |i| < & — >■ oo. (1.6) 

These heuristic considerations suggest that the optimal energy decay estimate in Theorem 1.1, modulo 
an e-loss, should be of the form 

E{u,t) 1 ' 2 <^0\\(u Q , Ul )\\ D{A) , e>0. (1.7) 

(Remark that taking t = ek, e -C 1 in (|1.6j) and letting k tend to infinity shows that (|1.7|) cannot hold 
if e < 0.) 

In section 3 we shall describe a class of examples of C°°-damping coefficients a for which we are 
able to improve the result of Theorem 1.1 and obtain an estimate of the form (|1.7|) . 



Theorem 1.2 Assume that a E C°°(f2) is a damping coefficient satisfying the assumptions of the 

beginning of section^ and in particular (|3.1j) . (|3.2[) . for some m > 4. Then we have for k > and 
e>0, 

^(w s *) 1/2 <^Pr-||(«o,wi)lb(>i*), t>2. (1.8) 
t^ £ 

In particular if (|3,2|) is satisfied for all m > 4, then for each k > and each e > there exists Ck, e > 
such that for each {uq,u{) G D{A k ) we have 

^(«,<) 1/2 <^fc=|||(«o,«i)||jj(>l*) J t> 2 - (1-9) 



At the present time, we do not know whether one can obtain an estimate of the form (|1.7|) for a 
general continuous, or even smooth, damping term, concentrated near the non-rectangular part of the 
domain. 

Remark. When this paper was complete, the authors learned of a very recent work by Kim Dang 
Phung jHj, where a polynomial decay estimate for the energy, with an unspecified decay rate, was 
obtained on a partially cubic domain, with the damping acting in a neighborhood of the boundary 
except between a pair of parallel square faces of the cube. An initial glance at [S] shows that the 
methods of the present paper are completely different. 

Acknowledgment. The second author is grateful to the National Science Foundation for partial 
support under the grant DMS-0304970. He would also like to thank Rowan Killip for a stimulating 
discussion. 



3 



2 Proof of Theorem 1.1 

Following [I] and [Rj, we shall use the stationary methods, and the main point will be to estimate the 
resolvent (A - .4)" 1 : H -► H for A G R, |A| > 1. 

The following result, closely related to Proposition 6.1 of [I], is our starting point. 

Proposition 2.1 Let R = [0, l] x x [0,7r] y C R 2 be a rectangle, and let us consider the stationary 
problem 

(-A + 2iaX -\ 2 )u = f + 8 x g, u\ dR = 0. (2.1) 
Here A € R, |A| > 1 and f,g£ L 2 (R). Assume that < a G C(R) is such that 

[0, 1]9ih inf a(x, y) does not vanish identically. (2-2) 

y6[0,7r] 

Then we have 

a(x,y) \u\ 2 dxdy) . (2-3) 



\u\\l 2(R) <0(l) (\\f\\h iR) + \\g\\h iR) +\ 2 j R < 



Proof: When establishing (|2,3j) . it suffices to do so when a = a{x) is a function of x only. Indeed, 
assume that (|2.3|) has already been established in this special case. If a G C(R) satisfies (|2.2j) . we can 
take < a\ G C([0, 1]), not identically zero and such that < a\{x) < a(x,y) when (x,y) £ R. Now 
if u G Hq (R) satisfies 

(-A + 2iAa - A 2 ) u = f + d x g, (2.4) 



then 

i2 



(-A + 2iAai - A 2 ) u = f + d x g + 2i\(a 1 - a)u. (2.5) 
Applying (|23|) to (|23|l . we get 



< O(l) ^|| / Ilia + ||s Ilia + A 2 ^ (a — ai) |n| 2 dxdy + A 2 J a± \u\ 2 dxdy^j , 



(2.6) 



and the bound ()2.3|) follows in the general case. In what follows, when proving H2.3|) . we shall therefore 
assume that < a = a(x) G C([0, 1]) is a function of x only, which does not vanish identically. 

When analyzing H2.1j) . we follow [1], [5] and separate variables. The Dirichlet realization of —d 2 
on L 2 ((0,7r)) has the eigenvalues k 2 , k = 1,2,... with the corresponding eigenfunctions efc(y) = 
\/2/ir sin ky, forming an orthonormal basis in L 2 ((0,7r)). Writing 

oo 

u{x,y) = ^2u k (x)e k {y), 

k=l 

and similarly for / and g, 

oo oo 

f(x,y) = ^2fk(x)e k (y), g(x,y) = ^ g k (x)e k (y), 

k=l k=l 

we see that u k , k = 1, 2, . . . satisfy 

{-d 2 x + 2ia(x)X + k 2 - X 2 )u k (x) = f k (x) + d x9k {x), u k (0) = u k (l) = 0. (2.7) 
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When analyzing (|2.7|l . we shall first consider the case when k G N is such that k > |A|. We then 
claim that 

IK Ilia < O(l) (|| Mlia + 11^11^), (2.8) 

and here we do not need the assumption that a does not vanish identically. 

Indeed, multiplying the first equation in (|2,7|) by u^, we obtain, integrating by parts and taking 
the real part, that 



\u' k {x)\ dx + {k 2 -\ 2 ) / \u k {x)\ 2 dx (2.9) 

Jo 

= Re / (fk(x)uk(x) - gk(x)d x Uk(x)) dx. 
Jo 

Here the right hand side does not exceed 

(II fk Hl2 + \\g k || L 2) (|| Wfc||z,2 + \\u' k \\ L2 ) < 2(11 f k \\ L 2 + HfffcH^a) \\u' k \\ L 2, 

where we also used the Poincare inequality 

IKIU 2 < ||u' fe || L 2, u k G £^((0,1)). 

Another application of the Poincare inequality combined with (|2.9jl then immediately gives (|2.8|) . 

We now come to discuss the second case when k < |A|. Let r G R be such that r 2 = A 2 — k 2 . Then 
we have 

(-<9 2 + 2ia(x)r - r 2 ) u k (x) = f k (x) + d x g k (x) + 2i(r - A)a(x)u fc (x), (2.10) 

u fc (0) =« fc (l) =0. 

Now the operator in the left hand side of (|2.1U|) is a one-dimensional stationary damped wave operator, 
for which it is well-known |B] that the resolvent bound 

Mr) := + 2iar - r 2 )' 1 = O (j^j ■ L 2 ((0, 1)) - L 2 ((0, 1)) (2.11) 

holds true, since < a £ C([0, 1]) is not identically zero. See also [7j. It follows that 

fl (r)=O(l):L 2 ((0,l))^F 1 ((0,l)), r G R, (2.12) 

and hence by duality, the same bound holds when Rq(t) is viewed as an operator from i? -1 ((0, 1)) to 
L 2 ((0, 1)). Applying these observations to ()2.1Uj) . we get 

\\u k \\ 2 L2 < 0(1) (|| A Hia + \\g k Ilia) + 0(1)A 2 || a^K ||i 2 . (2.13) 

Combining the bounds Q2.8|) and (|2. 13|) and summing with respect to fe, we get the proposition. □ 

Proposition 2.2 Let Q C R be a partially rectangular domain, £1 = RiJ W , and let < a G C(Q) be 
such that a > in W. Then if we put R(X) = (—A + 2zaA — A 2 )^ 1 , A G R, then we have 

R(\) = 0{\X\):L 2 (n)^L 2 (n), |A| > 1. (2.14) 
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Proof: We may assume that R = [0, 1} X x [0, ir] y , with the sides parallel to the x-axis contained in the 
boundary of fi. 

When / G L 2 (Q), let u £ Hq(Q) be the solution of 

(-A + 2mA - A 2 ) u = f, u\ 9n = 0, |A| > 1. (2.15) 

We now let < x ^ Co°((0, 1)), < x < 1> be a cut-off function such that x = 1 on [e, 1 — e]. Here 
e > small is to be chosen. Then it follows that x u i viewed as a function on R, vanishes on its 
boundary and satisfies in the interior 

(-A + 2mA - A 2 ) X u = xf+ [-A, xk (2.16) 
Here [— A,x]ti = x"( x ) u ~ 2dx(x'( x )u). An application of Proposition 2.1 gives therefore that 

||xu||| 2 <0(l)(||/|||2+ / \u\ 2 dxdy + \X\ 2 [ a(x,y)\u\ 2 dxdy] . (2.17) 

V Jw f JR J 

Here iv e is a neighborhood of the support of x'( x ), an d choosing e > small enough, we can achieve 
that u> e is so close to the vertical sides of dR, that it is contained in the set where a(x, y) is bounded 
away from 0. It follows that 



Xn||| 2 <0(1)||/||| 2 +0(1)|A| 2 f a(x,y)\u(x,y)\ 2 dxdy. 

Jn 



(2.18) 



Now observe that multiplying (|2.15|) by u, we obtain, integrating by parts and taking the imaginary 
part, that 

2A / a \u\ 2 dx dy = Im / fudxdy. (2-19) 
Jn Jn 

Hence 

|A| / a\u\ 2 dxdy <\\f\\ L 2\\u\\ L 2. (2.20) 
Jn 

and using this in (|2.18|) we see that 

\\xu\\ 2 L 2<0(l)\\f\\l 2 + 0(l)\X\\\f\\ L 2\\u\\ L2 . (2.21) 

It remains to estimate the L 2 -norm of (1 — x) u i an d to that end we notice that the support of this 
function is contained in the set where a is bounded away from zero. Therefore, another application of 
(EHUl) shows that 

||(1-XH|£ 2 <0(1)||/|U2|M| L3 , |A|>1. (2.22) 
Putting together the estimates (|2.21l) and (|2.22l) we get 

<0(l)ll/llL 2 +0(l)|A| 1/2 ||/||^ 2 ||n||^ 2 , (2.23) 

so that 

|M| L2 <0(1)|A|||/|U 2 . (2.24) 

This completes the proof. □ 

With Proposition 2.2 available, we are in the position to estimate the norm of the resolvent of 
the operator A in (jl.2|) on the real axis. In doing so, we notice that a simple computation, as in jS] 
and 0, shows that 
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Here, as in Proposition 2.2, we have written 12(A) = (—A + 2iXa — A 2 ) . To estimate the norm 
(A — .A) -1 as an operator on Tt = Hq x L 2 , we have to derive bounds on the operators 

12(A) : L 2 ^ H \ R(X)(2iaX - A 2 ) - 1 : H\ -► L 2 ', (2.26) 

and 

R(X){2ia - A) : H% -> Hq 1 . (2.27) 
Now Proposition 2.2 together with an integration by parts argument shows that 

12(A) = 0(X 2 ) : L 2 -» H \ (2.28) 

and hence by duality, 

12(A) = 0(A 2 ) : H' 1 -► L 2 . (2.29) 

When estimating the norm of the second operator in ()2.26|) . we use that 12(A) (2iaX — A 2 ) — 1 = 12(A) A, 
and hence, combining (|2.29|) together with the fact that the Laplacian A : Hq —* H _1 is continuous, 
we get 

R(X)(2iaX - A 2 ) - 1 = 12(A)A = 0{X 2 ) : Hq 1 -» L 2 . (2.30) 
It remains to estimate the norm in (|2.27|) . and to that end, we write 

12(A) (2m — A) = -r(l + 12(A) A). (2.31) 
A 

If / € Hi and u = 12(A) A/ G Hq 1 then 

(-A + 2mA-A 2 )u = A/er 1 , (2.32) 
and hence multiplying ()2.32j) by integrating by parts, and taking the real part, we get 

II "I lag ~ a2 IMIl2 < II A/Htf-ilMl^i < ||/||hi||u|| h i. (2.33) 

Therefore, 

\\u\\ 2 H , <0(1)(A 2 || U || 2 2 + (2.34) 

When estimating the L 2 -norm of u = R(X)Af, we use (|2.3U|) to conclude that 

\\u\\ L 2<0(X 2 )\\f\\ H i, (2.35) 

and combining this estimate with (|2.34|) we get 

IMIfli <0(|A| 3 )||/|| Hol . (2.36) 

An application of (|2.31|) then shows that 

R(X)(2ia - A) = 0(X 2 ) : Hq 1 -> Hq 1 . (2.37) 

Combining Proposition 2.2 together with (|2.28|) . ()2.3(Jj) . ()2.37|) . and the fact that A has no real 
eigenvalues, we get the basic bound 

(A - A)' 1 = + |A|) 2 ) : H -> W, A G R. (2.38) 
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We shall finally show how the bound (|2.38f) allows us to conclude the proof of Theorem 1.1. In 
doing so, we shall follow the argument of 8^ closely, adapting it to the present case. 
When k > 1 is an integer and x £ Ti, we write, as in [B] and for t > 0, 

e **(i -ur**= iriJ/ tx (T^ {x - A) ~ lxdx - (2 - 39) 

Here 7 = {AeC;A = r/ — i/2, rj £ R}. Furthermore, when X = X(t) = 71 (i/ logt) 1//2 , where 71 > 
is to be chosen, we use the same decomposition of ()2.39|) as in [Sj, 

e itA (l - iA)~ k x (2.40) 



1 1 r r itx 1 



e -, — -re 



2m J J R (1 - iX) k 



-(A-r) 2 /2( A _^)-l xdrdA 



••• + / / •••=: h+h, 

~{J\t\<X JjJ\t\>X 



where 



and 



1 1 r r 1 



C HX L c -(r-Xy/2 



2m yj2~TT J-y J\t\<x (l-iX) k 



{X-A)~ l xdrdX, (2.41) 



/ 2 = J_ * / f e «A 1 - (T -A)V2 (A _ A r l xdrd\. (2.42) 
2m^JJ lrl > x (l-i\) k v ; v 7 

Now combining 1)2.38 J) together with a perturbation argument, we see that the function A 1— > 
(A — ^4) _1 x is holomorphic in A with values in Ti, in the region below and including the curve 

A = r, + i% |r/| > ll U {A = t? + ieo, \v\ < 1} , (2-43) 



where eo > is small enough but fixed. Also, we see that along 7 eo , (A — A)~ 1 x is bounded in Ti by 
C(1)(1 + |t/|) 2 ||2;||^. Therefore, when estimating I\, we may write 

h = -?-^=[ [ e ux —^—re^-^/ 2 {X-A)- l xdrdX. (2.44) 
•Ziri^ J leo J\ T \<x (l-^A) fc V ; K ' 

The contribution to I\ coming from the part of 7 eo where \rj\ < 1 is easily seen to be bounded by 
(D(e- eot X(t))\\x\\ H , hence decaying exponentially, and therefore we may concentrate on the part of 
7 eo where \r)\ > 1. When estimating the contribution coming from this part, exactly as in 0, we 
treat separately the cases when 1 < \rj\ < 72 (i/ log i) 1 / 2 and |r/| > 72 (i/ logi) 1 / 2 . Here 72 > is to be 
chosen. The "H-norm of the integrand in I\ coming the part where \rj\ > 1 is bounded by 

0(1)ex K-fF)(rT^ e " ,M ' ,2/2||l|l «- ,2 - 45) 

Here |r| < 7i(t/logi) 1 / 2 . It follows that the contribution to I\ coming from the region where 1 < 
\v\ < 72 (^/ log t) 1 / 2 is controlled by 

0{1)t^ Wll^llwi ( 2 - 46 ) 
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which decays as any fixed inverse power of t, provided that 72 > is sufficiently small. Having fixed 
72 > small enough, as in 8], we choose 71 G (0,72), and observe that for |r| < 71 (i/ log i) 1 / 2 and 
> 72 (V logi) 1 / 2 , we have 

^(t-t^^ + t? 2 ), <5>0. (2.47) 

It follows then that the final contribution to I\ coming from the region where > 72 (t/ logi) 1 / 2 is 
obtained by integrating over this region 

and this clearly decays rapidly (and even exponentially) as t — > 00. We conclude that for each JVeN, 

h = N {t- N )\\x\\ n . (2.49) 
When estimating I2 in (|2.42j) . we continue to follow ^ and write 

I 2 = e itA J, (2.50) 

where 



J = _i__^/ / - \^ k e-^ 2 / 2 {\-A)- l xdTd\. (2.51) 



1 1 

27ri ^/2tt 7| r |>x J-y (! - ^) 

Since e**" 4 is uniformly bounded on H. for t > 0, it suffices to estimate J in the energy norm, as t — > 00 
Exactly as in [S], we then see that we have to estimate the integral 



oc 



S x _i +lle -i*/8(x)dn, (2.52) 
^0 2 

where 



S T (x) = J {l \ x)k ^ (X ~ T? (A - A)~ l x dX. (2.53) 

Arguing precisely as in jHj and making a contour deformation as in that paper, we then verify that 
with r = X - § + /ie- i7r / 8 , /x > 0, 

II^WII«<0(l)(«-*^' + IIT ^^ F )lb||„, 5>0. (2.54) 

Let us remark here that when deriving l|2.54|) . the bound (|1.3j) is important. 

It follows from (|2.5Uj) . 1)2.52(1 . and (|2.54l) that the W-norm of I2 does not exceed a constant times 
the norm of the vector x multiplied by 

Since 



/ 

Jo 



/„ Ou + l)( M + X)* V ^ fc 

we conclude, recalling the definition of X, that 

(logf) fc / 2+1 

I 2 = O k (iy-^j- 2 ||s|| w , t>2, fe = 2,3, ... (2.57) 

This completes the proof of Theorem 1.1 in the case when k > 1 is an integer. Using an interpo- 
lation argument as explained in [5], we get the result for a general > 0. 

Remark. By refining again the analysis above as in Section 4] one could probably avoid the 
logarithmic loss. 
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3 Improved decay estimates: an example 



In this section, we shall give an example of a class of damping coefficients vanishing away from a 
neighborhood of the non-rectangular part of the domain, for which the result of Theorem 1.1 can be 
improved, leading to an estimate of the type p. 7(1 . When doing so, to fix the ideas, we shall let = S 
be the Bunimovich stadium, defined in the example preceding the formulation of Theorem 1.1. 

When S = R U W, R = [0, l] x x [0,7r] v , we let < a e C°°(S) be such that a > in W. We 
assume that 

a-^O) = [5, (1 - S)] x x [0, n] y , < 6 < 1, (3.1) 

and that a is independent of y when x is close to 5 and 1 — 5. Furthermore, let us assume for simplicity 
that close to 5 (respectively 1 — 5), we have 

d m 

a(x) < (respectively > 0), (3.2) 



dx 

for some m > 4. We then immediately deduce the following result. 
Lemma 3.1 For any n < m there exists C n ^ m > such that 

a (n) (x)| <C n , m a(x)^. 
Proof: It suffices to consider the case x > (1 — 5). Using Taylor's formula, we get for x > (1 — 5), 

f-x I _ a \m-l f-x ( _ „\m—2 

a(x)= [X S) a^(s)ds, a'(x)= {X S) a^(s)ds. 
J is Jis (m-2)! 



As a consequence, 



\a'(x)\< ^ ^ aW(s)ds+ / ^ *> aW(s)ds, (3.3) 

J (1-5) (m-2)l J y (m-2)! 

with (1 — 5) < y < x to be chosen. To estimate the first integral in ((3.3|) we use the bound 

(x - s) m ~ 2 _ (m-l)(x- s) m ~ l < (m - 1) (x - s)" 1 " 1 



(m — 2)! x — s (m — 1)! x — y (m — 1! 

and for the second integral, we use that 

(x - s) m ~ 2 (x - y) m ~ 2 
(m - 2)! ~ (m-2)! 

We obtain 

|o'(x)| < ^-^a(x) + O m (l)(x - y) m ~\ 
(x - y) 

Choosing y so that a(x) = (x — y) m gives the lemma for n = 1. The general case is similar. □ 
Let us consider now the stationary problem 

(-A + 2iAa- \ 2 )u = f G L 2 (S), u\ ds = 0, A > 1. (3.4) 
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As in the proof of Proposition 2.2, we let < xx £ Co°((0> 1)) be a cut-off function so that xx vanishes 
for x close to or l.Then the function xx u vanishes on dR and satisfies in the interior of the rectangle, 



(-A - A 2 ) x\u = xf + Xx u - 2d x (x' x u) - 2i\a(x)x\u. 



(3.5) 



We choose xx = x(^ a ( x )) with x = for \x\ > 2, and x = 1 f° r l x l ^ 1> so that in the support of 
X^(ic) we have 

a(x)~ T , A>1. (3.6) 
A 

From Lemma 3.1 we get the following bounds on the derivative of XA; 

|xi| = |A X / (Aa(x))a'(x)|<0(l)A^, 



(3.7) 



and similarly 



(n) 



<C(1)A™, n < m. 



(3.8) 



We now come to estimate the L 2 -norms of the functions x'\ u i X\ u i an d \a(x)x\u, occurring in the 
right hand side of ()3.5|) . In doing so, we write 



X\u 



x'a l l 2 u 



and using ()3.6|) . I|3.7jl . we get 



and consequently 



Xx 



lV2 



< 0(1)A2- 



X>lli 2 (fl) < 0(l)AS + m || a}/ 2 u \\ L 2 {S) 



Estimating the L 2 -norm of x'\U in a similar way, we get 



X '{u\\ L 2 < 0{l)X^\\ a y 2 u\\ L 2 {s) . 



Finally, a similar argument shows that 

\\a\ X xu\\ L 2 {R) < 0(l)\V 2 \\aV 2 u\\ L 2 {S) . 



(3.9) 
(3.10) 
(3.11) 



At this point we can apply Proposition 6.1 of 4 to ()3.5|) directly, choosing to x there to be a 
neighborhood of the edges x = and x = 1 where a is bounded from below and consequently xx u 
vanishes. Using (|3.9|) . (|3.1(J|) . and ()3.11|) . we get, with u = lo x x [0, ir] y , 



\Xxu\\ 2 L2{R) <0(l)(\\f\\ 2 L2{R) 



+ A || a~'~u 



An application of (12.20(1 in (|3.12j) gives next 

\\XXU\\ 2 L 2 < 0(1) (\\ f \\ 2 L2 + \^\\ f \\ L 2\\ U \\ L 2 



(3.12) 



(3.13) 



It remains to control the L 2 -norm of (1 — Xx) u i an d when doing so we remark that the support of 
1 — xx is contained in the set where a > 1/A, A>1. Using (|2.20j) for the second time, we infer that 



(l-Xx)u\\ 2 L 2<0(l)\\f\\ L 2\\ U \\ L 2. 



(3.14) 
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Putting together Q3.13|) and (|3.14f) we get the estimate 

\\u\\ L 2 <0(l)(||/|| L2 + A^||/||^ 2 ||n||^ 2 ), (3.15) 

and finally, 

|| u \\ L 2 < 0(l)\\\^\\f\\ L i, AeR, |A|»1. (3.16) 

Repeating the arguments of section 2, with the bound (|3.16|) in place of Proposition 2.2, we get 
the bound 

(A - .A) -1 = 0(1) (1 + |A|) 1+ ™ :H^H, AeR, (3.17) 

and Theorem 1.2 follows. 
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